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We derive the steady state solution of the Fokker-Planck equation that describes the dynamics 
of the nondegenerate optical parametric oscillator in the truncated Wigner representation of the 
density operator. We assume that the pump mode is strongly damped, which permits its adiabatic 
elimination. When the elimination is correctly executed, the resulting stochastic equations contain 
multiplicative noise terms, and do not admit a potential solution. However, we develop an heuristic 
scheme leading to a satisfactory steady-state solution. This provides a clear view of the intracavity 
two-mode entangled state vaUd in all operating regimes of the OPO. A nongaussian distribution is 
obtained for the above threshold solution. 



I. INTRODUCTION 

Among the quasiprobability functions that represent 
the density operator of a quantum state, the Wigner dis- 
tribution has undoubtedly advantages over others, since 
in this phase space representation the quantum-classical 
correspondence is, in general, much more visible. At the 
same time these functions contain all information avail- 
able in the density operator. 

Unfortunately, just a few states have been repre- 
sented by this distribution due to the difficulty to solve 
Fokker-Planck equations for nonlinear systems. Indeed, 
the phase space description of quantum systems is well 
known for quadratic hamiltonians JLi 2\ , but very little 
is known for nonlinear systems, such as, for example, 
the single- mode degenerate parametric oscillator 3*, or 
the transverse multimode degenerate parametric oscilla- 
tor [3. 

Of special interest is the two-mode entangled state 
generated in the optical parametric oscillator (OPO) Q, 
used in many experiments of quantum information and 
recently shown to produce tripartite continuous variable 
entanglement [g . The Wigner distribution of this state 
was obtained exactly using the solution derived from the 
complex P-representation (after adiabatic elimination of 
the pump mode) Q, and tranformed into a Wigner dis- 
tribution. Since this analytical result appears as an in- 
finite series of gamma functions (hypergeometric series), 
some important physical aspects are hidden by its math- 
ematical complexity. For example, two-mode entangle- 
ment is not of easy identification from the usual partial 
transpose criterion [^, Q . The same happens to the gen- 
eralized P-representation for the intracavity parametric 
oscillator derived in Ref. [10] • In order to calculate the 
intensity correlations, the authors obtained an integral 
expression that depends on degenerate hypergeometric 
functions, and the outcome also depends on an infinite 
sum. 

In the present work, we derive the Wigner function of 
a two-mode quantum state of the electromagnetic field 
generated by an OPO using a truncated Wigner equa- 



tion. This truncated approximation is valid when the 
nonlinearity is small enough, which is true for most of 
the OPO experiments developed so far. A simple expres- 
sion is obtained for the Wigner distribution, which ren- 
ders clear the many quantum features of the intracavity 
OPO field such as squeezing and entanglement. More- 
over, this expression is valid in all regions of the OPO 
operation regimes, that is, below, and above threshold, 
where the validity of both the linear approximation, and 
the perturbation theory break down. This distribution 
also describes the OPO quantum fluctuations properly 
around the operation threshold, while the usual linear 
approach provides divergent results. In order to illustrate 
the reliability of the distribution we obtain, we compare 
moments calculated with this Wigner function with those 
obtained by numerical simulations using the positive P- 
representation of the density operator without any ap- 
proximations, that is, full quantum result for any control 
parameters of the OPO. 

The paper is organized as follows. Section lU presents 
the master equation describing the three-mode OPO dy- 
namics in the Wigner representation. The correct adi- 
abatic elimination of the pump mode is addressed in 
Sec. mil where we show that the standard procedure [il| 
leads to wrong results for some quadrature moments. 
After truncation and adiabatic elimination a two-mode 
Fokker-Planck is obtained which does not satisfy the con- 
ditions for a steady-state potential solution. However, 
we show that an approximately potential solution can be 
found which reproduces satisfactorily many quantum fea- 
tures of the OPO (Sec. lIV[) . Comparisons with the linear 
theory are exhibited in Sec. |Vl Section IVII contains the 
concluding remarks. 



II. THEORETICAL MODEL 

We present here a model of three quantized modes cou- 
pled by a nonlinear crystal inside a triply resonant Fabry- 
Perot cavity. The Heisenberg picture Hamiltonian that 



describes this open system is given by [l| 
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Here E represents the external coherent driving pump 
field at frequency loq. The operators 00,01 and 02 rep- 
resent the pump, signal and idler fields, respectively, 
satisfying the following frequency matching condition, 
loq = oji + LU2- The terms Ti represent damping reser- 
voir operators, and x is the nonlinear coupling constant 
due to the second order polarizability of the nonlinear 
crystal. 

The master equation for the reduced density operator, 
after the elimination of the heat bath by standard tech- 
niques 1], is given by 
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where 7^ is the corresponding mode damping rate. 

In order to treat the operators evolution, we now turn 
to the method of operator representation theory. These 
techniques can be used to transform the density ma- 
trix equation of motion into c-number Fokker-Planck or 
stochastic equations. The phase space Wigner equation 
for the nondegenerate parametric amplifier that corre- 
sponds to the master equation ^ is then 
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This is not a Fokker-Planck equation due to the third 
order derivative term, but in the case where x is small 
enough we can drop this term, as discussed in appendix 
A (see Ref. [12] for a recent review of the truncated 
Wigner approximation and its applications). The trun- 
cated equation so obtained is a genuine Fokker-Planck 
equation with a positive diffusion term, and we can eas- 
ily derive the corresponding stochastic differential equa- 
tions in the rotating frame {olj = Uj exp {—iujjt) with 
J =0,1, 2): 
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III. VALIDITY OF THE ADIABATIC 
APPROXIMATION 

It is possible to find a stationary solution of the 
truncated Fokker-Planck equation when we adiabatically 
eliminate the pump mode variables. This means that we 
are considering the relaxation rate of this mode much 
larger than those of the downconvertcd modes, that is 
7o ^ 7i , 72 • The stationary solution for the pumped 
mode is 

ao^ —[E- xoiia2 + Vlo^o{t)] , (5) 

70 

where the noise term is retained in the adiabatic elimi- 
nation in order to properly deal with the noise dynamics. 
We are then left with two nonlinear dynamical equations 
for the complex amplitudes of the down-converted fields, 
where the pump amplitude is replaced by expression ([5]) . 
We now define the following real quadrature variables: 
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so that the remaining dynamical equations can be cast 
in the compact form: 
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where X = [xi,yi,X2,y2]'^ is a column vector, and the 
drift vector is defined as 
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We have set 71 =72 =7, which is a reasonable physical 
assumption for most OPO experiments, and defined the 



normalized pump parameter fi — x-E'/(77o), as well as 
the nonlinear coupling g — x/iV^JJo)- B is a 4 x 6 
matrix defined as 
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and ^{t) is a six component column vector whose entries 
are uncorrelated real Gaussian noises associated with the 
noise terms for the three interacting fields. 

Note that the last two columns in B give rise to multi- 
plicative noise terms in the stochastic equation ([7]). It is 
important to remark that if these multiplicative noise 
terms are drop [ll| . some essential quantum features, 
particularly important in the above threshold regime, will 
be lost. 

In order to illustrate the importance of the multiplica- 
tive noise terms around and above threshold, and study 
the validity of both the truncation and adiabatic ap- 
proximations, we compare the numerical integration of 
the stochastic equations in the positive-P representation 
(without adiabatic elimination) with the stochastic equa- 
tions for the truncated Wigner representation in the adi- 
abatic approximation both with and without multiplica- 
tive noise terms. The results for the squeezed quadra- 
ture variances (defined in Sec. |V| as a function of the 
pump level is presented in Fig. [T] The positive-P calcu- 
lation clearly agrees with the truncated Wigner if the 
multiplicative noise terms are properly taken into ac- 
count. Although valid below threshold, neglection of the 
multiplicative noise terms clearly fails around and above 
threshold, where it predicts the unphysical 1131 1 vanishing 
oft 




Figure 1: Squeezed quadrature variance as a function of the 
pump level. We compare the results of the positive-P (solid 
black line) and truncated- Wigner stochastic equations with 
(solid red line) and without (dotted black line) multiplicative 
noise terms. The first two curves are almost indistinguishable; 
the third one decreases monotonically. We chose g^ = 0.01 
in all calculations and 70 = 10 7 in the positive-P simulation. 
All plotted quantities are dimensionless. 



IV. APPROXIMATE POTENTIAL SOLUTION 
OF THE FOKKER-PLANCK EQUATION 

The set of stochastic equations (O can be mapped onto 
a genuine Fokker-Planck equation: 



dWjX.) 
dt 



d i_d d_ 



W{X.) (10) 



where the diffusion matrix is defined as D = BB"^, and 
summation over repeated indices is assumed throughout 
the text [IJ]. 

Neglection of the multiplicative noise terms in the adi- 
abatic elimination allows for a steady state solution in 
the potential form W{X) = TVexp (- / Z., dXi) with Z^ 
given by 
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as obtained in Ref. [TlJ. As we have seen, this distribu- 
tion gives incorrect squeezing results for above threshold 
operation. On the other side, if multiplicative noise terms 
are kept, a potential solution is no longer possible 

This difficulty can be circumvented if one replaces the 
diffusion matrix by its mean value, as explained in Ap- 
pendix B, so that a potential solution for W(K) becomes 
possible. This solution provides a quite simple form for 
the Wigner distribution which allows for the calculation 
of several statistical properties. We remark that this 
approximate potential solution is achievable only when 
71 = 72 = 7 as we have already assumed |19j . In this 
case, the steady state Wigner distribution for signal and 
idler reads: 
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where Af is the normalization constant and unimportant 
terms 0(.g*) were neglected. We also defined 
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As a concrete physical example, expression ([T^ can 
describe the joint Wigner distribution for two polariza- 
tion modes of a frequency degenerate type II OPO, where 
signal and idler are distinguished by their polarization 
states. All statistical properties, including experimen- 
tally accessible quantities like quadrature noise and corre- 
lations in the stationary state may be calculated with this 
distribution in any operation regime of the OPO within 
the validity of the adiabatic elimination. It is interest- 
ing to note that the distribution given by expression (J12p 
is single peaked below threshold operation (fi < 1) and 
becomes double peaked above threshold. In Fig. 2, it is 
possible to visualize the conditional Wigner distribution 





Figure 2: Conditional Wigner distribution for y\ = and 
j/2 = 0, a) below (/i = 0.5); b) at {fi — 1.0); and c) above 
threshold {fj, = 1.5). In all cases g^ = 0.01. All plotted 
quantities are dimensionless. 



for fixed values of j/i and 2/2, for an OPO operating be- 
low, at, and above tlireshold. It is easy to identify the 
squeezed and unsqueezed quadratures in all regimes. 

For the nonlinear media employed in most OPO de- 
vices, the coupling parameter g <^ 1. Therefore, the 
first two terms in the exponent of the Wigner distribu- 
tion governs the OPO behaviour below threshold. They 
are enough to predict the usual nonclassical features. In 
fact, it is evident now that the usual linearized approach 
to quantum noise in below-threshold OPOs is equivalent 
to completely neglect the g'^ term in the distribution. 
However, it is important to notice that for /U > 1, the 
Wigner distribution becomes divergent if we neglect the 
g^ term. This means that this term plays a crucial role 



for a complete description of the OPO behaviour at and 
above threshold. 

Additional insight is provided by a careful investiga- 
tion of the marginal distribution obtained by integrating 
W{xi,yi,X2,y2) with respect to one mode variables. Of 
course, due to the symmetry of W, the form of the re- 
sulting marginal distribution is independent of the mode 
being traced out. The marginal distribution for mode 2 
becomes: 



W{x2,y2) = 



2nJ^fi 
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It gives the statistical properties of isolated measure- 
ments on mode 2. The variances given by this marginal 
distribution are larger than those of the vacuum state. 
This excess noise can be seen as a consequence of infor- 
mation loss when one looks only at part of the whole sys- 
tem. Fig. [3] presents the marginal distribution below, at 
and above threshold. Above threshold the marginal dis- 
tribution is clearly peaked out of the phase space origin, 
what is a consequence of the macroscopic amplification 
of the down-converted fields. However, the distribution 
presents radial symmetry indicating complete phase un- 
certainty in each individual mode. 



V. COMPARISON WITH THE LINEARIZED 
QUANTUM APPROACH 

Almost the totality of theoretical works devoted to the 
analysis of quantum noise in optical systems rely on lin- 
earization of the small quantum fluctuations around the 
macroscopic steady state mean values. It is worth to no- 
tice that this procedure has limited validity, especially 
around the threshold critical point where quantum fluc- 
tuations may become comparable to the mean values. In 
order to evidence the breakdown of the linearized theory, 
we shall compare its results with those obtained from the 
solution of the Fokker-Planck equation (fTO|) . 

The usual approach given to quantum noise in the 
literature is obtained from first order stochastic equa- 
tions of motion [l^ [I^ . These equations are often used 
to predict squeezing in a linearized fluctuation analy- 
sis. They are non-classical in the sense that they can 
describe states without a positive Glauber-Sudarshan P- 
distribution, but correspond to a Gaussian Wigner dis- 
tribution. 

We now find it useful to introduce combined field 
quadratures, as in two-mode approaches used previ- 
ously [l6| . These combined quadratures are the Einstein- 
Podolsky-Rosen (EPR) variables used to characterize 
continuous variable entanglement between the down con- 
verted fields. They are defined as 



a;± 



X\_ ± X2 



y± 



yi±y2 
V2 ■ 



(15) 






Figure 3: Marginal distributions: a) below threshold, fi = 0.8; 
b) at threshold, fi — 1.0; c) above threshold, ^ = 1.2. In all 
cases g^ = 0.01. All plotted quantities are dimensionless. 



These quantities correspond to the squeezed and anti- 
squeezed combined quadratures obtained in the lin- 
earized theory. From the first order stochastic equations 
one can easily obtain the steady state variances of the 
EPR variables fHfli: 
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for below threshold operation, and 
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for above threshold operation. 

Now let us briefly discuss the predictions of the lin- 
earized approach and its validity. The quadratures X- 
and J/+ exhibit the expected squeezing, going from the 
vacuum fluctuations for zero pump until 50% intracavity 
squeezing (which corresponds to perfect squeezing out- 
side the cavity) at and above the oscillation threshold. 
The unsqueezed quadratures x+ and y_ present diver- 
gent behaviour around threshold, which is certainly not 
physical. In fact, as we shall see next, the steady state 
solution of the Fokker-Planck equation (fTO|) gives a well 
behaved Wigner distribution which does not display any 
divergences. This Wigner distribution will also put limits 
on the amount of squeezing attainable. 

In order to investigate the two-mode entanglement di- 
rectly from the Wigner distribution, we now write it in 
terms of the EPR variables. Below threshold, where the 
linearization of the equation of motion is valid, we can 
neglect the g^ term, and the distribution can be approx- 
imated as 



WL(x+,y+,X-,y-) ^ Mexp 
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(18) 
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With this expression we can easily calculate any mo- 
ment of the distribution. For instance, we easily reobtain 
the results of Eqs. (fT6|) for the intracavity noise squeezing 
in the combined quadratures. As we mentioned before, 
at threshold we have perfect external squeezing while the 
unsqueezed combined quadratures blow up showing the 
failure of linearization [1^. Moreover, it is easy to char- 
acterize this state as entangled by using the Duan-Simon 
criterion [^, Q ■ Above threshold we can also calculate the 
moments and we find that this criterion shows the suffi- 
cient condition to characterize this state as entangled, in 
spite of not being a Gaussian state in that regime. 

It is also interesting to compare the results obtained for 
(x^) with the linearized theory (both below and above 
threshold) with those obtained from the Wigner func- 
tion given by Eq. dH]). This is done in Fig. (JH). While 
the linearized approach gives a divergent behavior at the 
oscillation threshold, the results obtained from the trun- 
cated Wigner function gives a more realistic smooth be- 
havior in a gree ment with the full quantum solution of 
Refs. [3, [lO, [201 • The results for the squeezed quadrature 
(x?_) are shown in Fig. [5] where, again, the results given 
by the truncated Wigner distribution agree with those 
obtained from the full quantum solution. Below thresh- 
old, the linearized approach also agrees with the full 
quantum theory and the truncated Wigner approach, but 
an important deviation can be observed above thresh- 
old where the linearized theory gives (x?.) = 0.5 for any 
pump power. The same results are obtained for the other 
squeezed variable (yi^), so that the Duan-Simon sepa- 
rability criterion [8|, |9| shows that the two-mode quan- 
tum state is entangled. However, note that both the full 




Figure 4: Comparison between the values of the anti-squeezed 
variable {x\) obtained from the linearized theory (dashed 
line) with those obtained from the Wigner function of ref- 
erence [31 (solid line), and the one given by Eq. (I12p (crosses). 
In all cases we have set g^ = 0.01. All plotted quantities are 
dime""'""'™" 
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Figure 5: Comparison between the values of the squeezed 
variable (a;?.) obtained from the linearized theory (solid line) 
with those obtained from the Wigner function of Ref. J,] (open 
circles), and the one given by Eq. (|12|l (rhombs). In all cases 



we set g = 0.01. All plotted quantities are dimensionless. 



VI. CONCLUSION 

By truncating the evolution equation for the Wigner 
function of an optical parametric oscillator, followed by 
the proper adiabatic elimination of the pump mode, we 
obtained a Fokker-Planck equation which describes sat- 
isfactorily the properties of the steady state in any oper- 
ation regime of the OPO within the validity of the adi- 
abatic elimination. This was checked by direct numeri- 
cal simulation of the corresponding stochastic equations. 
Then we showed that the Fokker-Planck admits a simple 
"quasi-potential" solution. For the unsqueezed quadra- 
tures, this solution does not present the unphysical diver- 
gent behavior of the linearized theory on the oscillation 
threshold. Moreover, it provides good agreement with 
exact quasiprobabilities already available in the litera- 
ture, both for the squeezed and unsqueezed quadratures. 
While these exact solutions are given as infinite series, 
the potential solution of the truncated Wigner distribu- 
tion presents a simple form allowing for an easier visual- 
ization of the phase space distribution. 
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APPENDIX A 

In order to clarify the hierarchy of the terms in the 
dynamical equation for the Wigner distribution, let us 
define the following dimensionless parameters: 

• Dimensionless coupling: 



quantum approach and the truncated Wigner distribu- 
tion predict a less pronounced violation of the criterion 
above threshold. 

In Refs. [lO, 1201 Fokker-Planck equations were derived 
for the complex-P [10] and positive-P [201 representa- 
tions, and stationary solutions were obtained for the 
corresponding distributions. Since these Fokker-Planck 
equations were derived without any truncation, we can 
consider the stationary distributions so obtained as ex- 
act. In Ref. [7[ an exact stationary Wigner function was 
mapped from the complex-P distribution of Ref. [iO| and 
expressed in terms of Bessel functions. It is important to 
remark that the simple Wigner function presented here 
provides a very good approximation for those exact so- 
lutions, as can be seen in Figs. [3| and [SJ even for a rather 
large value for the nonlinear coupling g^. 
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Rescaled field amplitudes: 
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with Wo = wi + aj2 (phase match) and Ec = 77o/x 
(threshold pump amplitude). 

Perturbative approaches to the quantum stochastic 
equations of the parametric oscillator has been widely 
discussed in Refs. [l^- In these approaches, the pertur- 
bation parameter is the dimensionless coupling g, usually 
assumed to be small as is the case in most experimental 
situations. The quantum fluctuations are then expanded 
in powers of g in order to provide progressive corrections 
to linearized theories (the zero order terms correspond to 
the macroscopic mean field equations). 

Therefore, let us rewrite the phase space Wigner equa- 
tion for the nondegenerate parametric amplifier in terms 
of the rescaled amplitudes and dimensionless parameters. 
The equation then becomes: 

dW ( d 
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This rescaled form, evidences the hierarchy of the 
different terms in the dynamical Wigner equation and 
makes it clear that the third order derivative terms are 
negligible for small coupling, since it is proportional to g*. 
It is interesting to make a close inspection in the physi- 
cal meaning of this hierarchy. For small enough coupling 
{g ^ 1), if we neglect both the g^ and 5^ terms, we are 
left with a Liouvillian evolution which is deterministic in 
essence, any indeterminacy comes from the initial condi- 
tions. In this case, for example, initially delta distributed 
amplitudes evolve as such for all times: 

W{t) = CS{f3o~Mt))S{f3i~Mt))S{f32-f32{t)). (27) 

This corresponds to a classical deterministic evolution of 
the system. 

As g is continuously increased, the next approximation 
corresponds to neglecting only the g^ term. In this case, 
Eq. ([26)) assumes the form of a standard Fokker-Planck 
equation, leading to the usual diffusive evolution with a 
possible stationary state (under self pulsing conditions, 
for instance, there is no stationary solution). 



Finally, there are effects of the g'^ term which are es- 
sentially quantum mechanical, but have received little 
attention in the literature. However, these effects require 
a coupling constant consideraby larger than what is ex- 
perimentally available in parametric systems. 



APPENDIX B 

Here we clarify the steady state solution derivation of 
the effective Fokker-Planck equation for the two down- 
converted modes of the field, after adiabatic elimination 
of the pump. We begin with the set of stochastic dif- 
ferential equations in the rotating frame, associated with 
the truncated Fokker-Planck equation for three modes 
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-7oao +E - xa-ia2 + \/7oCo(0 

-71^1 + xao"2 + \/7rCi(0 
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The adiabatic elimination of the pump mode requires a 
special care. We need to keep the fluctuations in the 
steady state beacause they provide important additional 
coupling between the down-converted modes. The sta- 
tionary solution of the pumped mode will be taken as 

ao = — (-B-xaia2 + 77^.^0 (i)) (29) 

70 

We now substitute this expression into the other two am- 
plitude mode equations to have the effective dynamical 
nonlinear equations, 

dai 1 

—— = -7iao + Xa2— (^-X"iQ:2 + V7o^o) + VTiti 
at 7o 

da2 * 1 / rn , — ^ \ , — ^ 

—- = -72ao + Xai— (^-X"ia2 + V7o^o) + V7242 
at 70 

(30) 

In the vacuum state (aja*) = 1/2, and {ii[t)*i][t')) = 
SijS{t — t'). It is clear from the above equations that the 
pump noise enhances the coupling between the ampli- 
tudes (quadratures) of the down-converted modes. 
In terms of real-quadrature variables. 



xi ^ ai+ a^ yi - 

X2 = a2 + Q*2 J/2 



-i (a.\ -\- aj) 
-i (a2 -I- 0*2) 



(31) 



the stochastic equations pop can be written in the form 

(32) 
where X is a vector that has the four quadratures as 



dX . „ , ^ 
-^A + B,(.) 



components, and the drift vector is defined as 

-xi + ^X2 - ^xi {x\ + y|) 



A = 7 



-y\ 


- ^J■y2 


X2 


+ ^J.Xl 


-y2 


-A^J/i 



-yi (a 



j/i) 



- ^X2 {xl + yl) 

- ^y2 {xl + yl) 



(33) 



if the potential condition, diZj = djZi , is satisfied. 
In the present case the Z vector is written as 



Zi = 



1 



l + 2^(xl+yf+xl+y^) 



-(1 + g^)xi + HX2 - —xx {x\ + y^) 



where we have defined the adimensional pump parameter 

= . B is a 



Z. 



/i — ^ — , and the nonlinear coupling g = 
4x6 matrix defined as 



B = v^ 



^10 ^x, ^y, ^^ 

10 -^V2 ^X2 



10 ^xi ^yi 



2 ^ V2 
9-5,, JL. 



(34) 



yOOO 1 -^yi ^xiy 



and r\(t) is a vector with six uncorrclated real gaussian 
noise components. 

Associated with this set of stochastic equation there is 
a Fokker-Planck equation that can be written as 



ap(x) 

dt 



d 

ax,. 



2 dX, 9X, '^ 



P(X) (35) 



where the diffusion matrix is defined as D = BB and 
has the following expression 



D 



Its inverse reads 



D 




(36) 



7 (ab - 



(P) 



b 



— c 
-d 



where we have defined 
a = 1-F 
b = 14 






yl) 



yl) 




(37) 



(a;iX2 +yiy2) 



l + ^ixf +yl+xl+yl) 



Z. = 



-(1 + 9^)yi + MJ/2 - yyi (a;? + yl) 
1 



^4 = 



i + ^(a;?+y? + a:i + yi) 

-(1 + .9^)3:2 + iixi - ■—X2 {xl + yl) 

1 

l + al(xl+yl+xl+yl) 

-(1 + 9^}y2 + Myi - yJ'2 (2:2 + 2/2) 



(40) 



Unfortunately, a potential solution cannot be achieved 
with the present form of matrix D. Of course, a steady 
state solution of the Fokker-Planck equation still exists, 
since the detailed balance is established, but the solution 
will not be in a potential form. 

For a sufficiently peaked distribution, we may approx- 
imate the diffusion matrix D by its average. The off 
diagonal terms have zero mean value due to phase space 
symmetry. Indeed, it is easy to verify that the dy- 
namical equations remain unchanged under the trans- 
formation a'l — ai exp(— 10) and a'2 = a2exp{iO). 
which corresponds to a definition of rotated quadra- 
tures Xi,y'i,x'2,y'2- This symmetry of the dynamical 
equations imply in VK(xi, yi, X2, j/2) = W{x[,y[,x'2,y2) 
For a particular choice of = 7r/2, this corresponds 



2/1: 



y'l 



-Xl, x'2 = — 2/2, and y'2 = X2 so that 
W{xi,yi,X2,y2) = W{yi,-xi,~y2,X2). This parity 
property of the Wigner function forces the off diagonal 
terms c and d to be zero. Note that the diagonal terms of 
the diffusion matrix are related to the individual intensi- 
ties of the down-converted fields and can be replaced by 
the mean values calculated from the deterministic classi- 
cal equations. This substantially simplifies the expression 
for vector Z 



d 



y {xiy2 - yiX2) 



(38) 



and ab- c^ - d'^ ^1 + S^ {xj+yf+ xl+ yl) 

The Fokker-Planck equation admits a potential steady 
state solution in the form -P(X) = iVexp — J^^ Z^X^, 
where Z^ given by 



E^i^ 



2Ak-J2 



ax, '^ 



(39) 



Zi 

Z2 
Z3 

Za 



1 



s(^) 

1 
s{[i) 

1 

s(^) 

1 

s(^) 



-Xl + [1X2 



-y\ + ms2 



-77Xl{^ 






y ( 2 

-X2 + [IXi - —X2 [X-i 



-2/2 + M2/1 - y 2/2 {xl 



yl) 



yl) 



vl) 



yl) 



(41) 



where s{fj,) = 1 + g"^ 



yf)/2 is given by Eq. ([TS 



This vector represents a conservative field that gener- 



ates the approximate potential solution of the effective 
Fokker-Planck equation. 
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